Structure of quantum disordered wave functions: weak localization, far tails, and 

mesoscopic transport 
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We report on the comprehensive numerical study of the fluctuation and correlation properties of 
wave functions in three-dimensional mesoscopic diffusive conductors. Several large sets of nanoscale 
samples with finite metallic conductance, modeled by an Anderson model with different strengths 
of diagonal box disorder, have been generated in order to investigate both small and large devia- 
tions (as well as the connection between them) of the distribution function of eigenstate amplitudes 
from the universal prediction of random matrix theory. We find that small, weak localization-type, 
deviations contain both diffusive contributions (determined by the bulk and boundary conditions 
dependent terms) and ballistic ones which are generated by electron dynamics below the length 
scale set by the mean free path £. By relating the extracted parameters of the functional form of 
nonperturbative deviations ( "far tails" ) to the exactly calculated transport properties of mesoscopic 
conductors, we compare our findings based on the full solution of the Schrodinger equation to difi'er- 
ent approximative analytical treatments. We find that statistics in the far tail can be explained by 
the exp-log-cube asymptotics (convincingly refuting the log-normal alternative), but with parame- 
ters whose dependence on £ is linear and, therefore, expected to be dominated by ballistic effects. It 
is demonstrated that both small deviations and far tails depend explicitly on the sample size — the 
remaining puzzle then is the evolution of the far tail parameters with the size of the conductor since 
short-scale physics is supposedly insensitive to the sample boundaries. 



PACS numbers; 73.23.-b, 72.15.Rn, 05.45.Mt, 05.40.-a 



I. INTRODUCTION 

Quantum coherence, its nonlocal features, and ran- 
domness of microscopic details can cause large fluctu- 
ations of physical quantities in disordered mesoscopic 
systems. The paradigmatic case is that of conductance 
fluctuations which has given impetus for the whole field 
of mesoscopic physics^'^ by pointing out at unexpected 
features of such fluctuations.'^ Contrary to the intu- 
ition developed from thermal fluctuations, and their self- 
averaging properties in the statistical physics of macro- 
scopic systems, the average value and variance are not 
enough to characterize the broad distribution functions of 
various mesoscopic quantities,"* even well into the metal- 
lic regime 3> 1 (5 = G/Gq being the dimension- 
less zero-temperature conductance, in units of conduc- 
tance quantum Gq — 2e^ /h). The fluctuations increase, 
broadening the distributions, as disorder is increased 
eventually driving a system through the localization- 
delocalization (LD) transition^ at g ~ 1. Thus, the meso- 
scopic program was born where full distribution functions 
of relevant quantities in open (e.g., conductance, local 
density of state, current relaxation times, etc.) or closed 
(e.g., eigenfunction amplitudes, polarizability, level cur- 
vatures, etc.) samples are to be studied. -^'^ Especially 
interesting are the large deviations of their asymptotic 
tails from the ubiquitous Gaussian distributions (which 
are expected only in the limit g — s- cx)). 

Recently, the study of fluctuations and correlations of 
eigenfunction amplitudes has been initiated.^ The quan- 
tum coherence induces long-range spatial correlations 



(due to massless modes, like diffusons and Cooperons) in 
the local density of states and eigenfunction amplitudes, 
which in turn lead to strong mesoscopic fluctuations of 
global quantities like conductance. Small deviations of 
eigenstate statistics from the universal predictions (ap- 
plicable in the limit g — ^ cxj) of random matrix theory 
(RMT) are well understood through perturbative correc- 
tions ^ 0{g^^) of the weak localization (WL) type,^ but 
the physical origin of large deviations in the far asymp- 
totic tail of the distribution function is much more con- 
troversial.®"^* Not only that there are different analytical 
predictions for the far tail asymptotics (which in fact do 
not explain all details of the tails found in numerical sim- 
ulations*^'*'^), but there is also an issue** of the relevant 
physics which is responsible for large wave function am- 
plitudes (quantum vs semiclassical) and a closely related 
question on the limitations of usually employed field- 
theoretical approaches*^ to study the disordered electron 
problem. Also, the parameters of the WL correction to 
the RMT framework cannot be explained (e.g., in the 
Anderson model*^) solely by the standard universal (in- 
dependent of the details of disorder) quantities extracted 
from the semiclassical diffusive dynamics.^ Instead, care- 
ful examination of ballistic effects, generated by the prop- 
erties of particle dynamics on the length scale below the 
mean free path is required.*^ Moreover, it is possi- 
ble that some types of disorder could generate apprecia- 
ble higher order terms, (characterizing non-Gaussian fea- 
tures of random potential*^) in this perturbative expan- 
sion in 1/(7, and thereby change the functional form of the 
perturbative correction as well. Thus, a detailed study of 
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deviations from the RMT statistics in the paradigmatic 
case of a quantum particle in a random potential offers 
a possibility to unravel underlying correlations in a con- 
trolled fashion, which paves the way for understanding 
plethora of diverse problems (including those outside of 
physics^^) where matrices containing random elements 
and their eigenstates are encountered. 

In the course of exploration of mesoscopic fluctuations, 
the so-called prelocalizcd states have been unearthed as 
the microscopic origin of asymptotic tails of various quan- 
tum distribution functions of thermodynamic and kinetic 
quantities. While typical wave fimction is spread 
uniformly throughout a metallic sample of volume L"^ 
with average amplitude L^'^^^ (up to inevitable Gaus- 
sian fluctuations), the prelocalized state in 3D exhibits 
much larger local amplitude splashes (on the top of the 
homogeneous backgromid |'I'(r)p ~ L~'^) at some points 
r within the sample. ^'^^ To obtain the far tail of such 
distribution "experimentally" (e.g., in microwave cavi- 
ties of Rcf. 21 or by numerical simulations'^'^) in "real- 
istic" metallic systems, one has to search for extremely 
rare disorder configurations where quantum interference 
effects are able to generate highly unusual eigenfunctions. 

What is the relevance of prelocalized states for trans- 
port experiments? Most of phenomena in good metallic 
disordered conductors are semiclassical in nature. This 
means that disorder-averaged properties, like conduc- 
tance measured in experiments or calculated in (quan- 
tum transport) theory,^^ are determined by the usual ex- 
tended states of uniform amplitude, formed in the typ- 
ical fluctuations of the random potential. However, re- 
cent experiments on quantum dots (nanofabricated sam- 
ples with well-resolved electron energy levels) show that 
some transport properties, like fluctuations of the tun- 
neling conductance, can depend sensitively on the local 
features of wave functions which couple the dot to exter- 
nal leads. ^'^ Also, to understand the excitation and ad- 
dition spectra of quantum dots one has to deal with the 
statistics of Coulomb interaction matrix elements, which 
are influenced by the eigenfunction amplitude fluctua- 
tions.^^ By exploiting the correspondence between the 
Schrodinger and Maxwell equations in microwave cavi- 
ties, it has become possible to probe directly the micro- 
scopic structure of wave functions in quantum disordered 
or quantum chaotic systems. 

Here we undertake a comprehensive search, through 
numerical simulations, for special disorder configurations 
in order to investigate functional dependence of the pa- 
rameters determining eigenstate statistics on the disorder 
strength or sample size. This is not just a 'brute force' 
study culminating in a fitting procedure of the observed 
distribution functions, but more importantly, an attempt 
to quantify those effects which can lead to substantial 
deviations from the RMT, departing even from the stan- 
dard semiclassical corrections to it. Namely, our results 
stem from the exact solutions of the Schrodinger equa- 
tion for a particle in a random potential, and therefore 
provide a reference point for the analytical approaches 



which usually integrate out some degrees of freedom by 
focusing on the "low energy" sector of a full theory.^ For 
this purpose, we also compute exactly the transport prop- 
erties of our finite-size samples, and relate them to the 
parameters extracted from the fits of analytical formulas 
to the perturbative and far tail interval of the eigenfunc- 
tion amplitudes. Mesoscopic physics intrinsically deals 
with finite-size phase-coherent samples, and has led to 
efficient use of different transport formalisms. Thus, we 
exploit the fact that transport properties of a specific 
sample (simulated here as nanoscale single band conduc- 
tors) can be measured exactly on a computer. Although 
our focus is primarily on the peculiar states exhibiting 
the largest amplitudes (which generate far tails of the 
statistics of eigenfunction amplitudes, as well as other 
mesoscopic quantities), it becomes necessary to investi- 
gate thoroughly the region of small eigenfunction am- 
plitudes because of the possibility that the same semi- 
classical quantities (like classical diffusion propagator, 
which we evaluate explicitly for the samples with specific 
boundary conditions) might govern both portions of the 
distribution function.'' 

We have investigated five different ensembles^^ of 
mesoscopic samples, each containing 30000 weakly dis- 
ordered three-dimensional (3D) metallic conductors. 
Finite-size samples are modeled by a tight-binding 
Hamiltonian (TEH) 

.H' = ^em|m)(m| -f ^ imn|m)(n|, (1) 

m (m,n> 

with nearest-neighbor hopping tmn = 1 (unit of energy) 
between s-orbitals (r|m) = tp{r — m) located on sites 
m of a simple cubic lattice of size L = 12a to L = 20a 
(a is the lattice spacing). Periodic boundary conditions 
are chosen in all directions. The disorder is simulated by 
taking potential energy to be a uniformly distributed 
random variable, —W/2 < £m < W/2, which is the 
standard (Anderson) model in the localization theory.^ 
The impurity configurations vary from sample to sample 
for a given disorder strength W G {3.5,4.0,4.5,5.0,5.5} 
which is chosen to ensure that ensemble-averaged trans- 
port quantities characterize metallic {g ^ 1) diffusive 
{£ <^ L) transport regime at half filling (i.e., at Fermi 
energy Ep = 0). For a weak kpt ^ 1 (fc^? is the Fermi 
wave vector) disorder, disorder-averaged transport prop- 
erties are semiclassical, i.e., well-described by the Bloch- 
Boltzmann formalism. Going beyond ~ 6 (but below 
critical Wc — 16.5 for the localization of the whole band) 
would still give metallic conductance g ^ 1 (for large 
enough lattice), but the semiclassical concepts appear- 
ing in analytical predictions for the far tail,'' like £, loose 
their meaning. The disorder strengths below = 3 
are excluded only because of requiring too large lattices 
to avoid quasiballistic transport. 

Statistical properties of eigenstates in a closed sample 
are described by a disorder-averaged distribution'°'^^ of 
eigenfunction "intensities" |'J'a(r)|^ 
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FIG. 1. Statistics of eigenfunction "intensities" |^'B(m)|^ 
in the band center E € (—0.15, 0.15) of an Anderson model on 
a simple cubic lattice 12'*. Each curve is obtained by examin- 
ing about two million exact eigenstates of the (time-reversal 
and spin-rotation invariant) Hamiltonians whose on-site po- 
tential random variable uniformly distributed over the inter- 
val [-W/2, W/2]. The disorder strengths W for the five dif- 
ferent sets of 30000 Hamiltonians arc chosen to ensure the dif- 
fusive (L S> i) and semiclassical {kpi > 1) transport regime 
in the conductors they model (relevant transport quantities 
are listed in Table I.) The random matrix theory prediction 
for the limit p — » oo is Porter-Thomas distribution, plotted 
here as a reference. 

= ~ |*a(r)| V)5(£; - E^)j , (2) 

on N discrete points r inside a sample of volume V. Here 
p{E) = d{E — Ea)) is the mean level density at en- 
ergy E, and (. . .) denotes disorder-averaging. Normaliza- 
tion of eigenstates gives t = J dtt f{t) = 1. We evaluate 
this function for eigenstates in the band center E = 0, 
which are obtained by exact numerical diagonalization 
of the Hamiltonian (1) [note that energy is a parame- 
ter in f{t)]. The distribution functions f{t) for the five 
sets of conductors modeled on the lattice 12^ is shown 
in Fig. 1. The explicit dependence of f{t) on the sam- 
ple size is demonstrated in Fig. 2 where W = 5 case 
is studied also on the 16"^ and 20^^ lattices. By searching 
through many configurations of the random potential one 
can find the rare ones which are responsible for the ap- 
pearance of states with the highest possible amplitude 
splashes. We smooth out the data by additional aver- 
aging over a small energy interval (which taken alone, 
or combined with only a small number of disorder real- 
izations, is not enough to study the prelocahzed states), 
without introducing any artifact in the computed distri- 
bution functions. This finally brings the number of an- 
alyzed eigenstates to about 2 • 10^ for each curve plotted 
in Fig. 1. 

The rest of the paper is organized as follows. In Sec. II 

we provide a short survey of the main analytical results 
for the eigenstate statistics in 3D, which are expected 
to be relevant for our observations. Then in Sec. Ill, a 




FIG. 2. Comparison of the W = 5 eigenstate statistics f{t) 
from Fig. 1 for a lattice 12'' and the statistics of band center 
eigenstates generated in the ensembles of 2000 or 500 conduc- 
tors modeled on larger lattices 16^ (with ~ 150 eigcnfunctions 
picked in a small energy interval AE = 0.3 around E = in 
each sample) or 20^ (with ~ 300 eigcnfunctions investigated 
in each sample), respectively. 

detailed comparison between these predictions, for both 
small and large deviations of distribution functions from 
RMT, and our results is undertaken using the calculated 
transport properties of the samples. Finally, we conclude 
in Sec. IV by looking beyond the raw numbers and point- 
ing out at open questions. 

II. STATISTICAL APPROACHES TO 
NONINTEGRABLE QUANTUM SYSTEMS 

Quantum dynamics of a non-interacting particle in 
random potential (e.g., generated by quenched impuri- 
ties) has a long history of being a standard playground 
for the development of ideas of Anderson localization^ 
and mesoscopic physics.^ The classical counterpart of 
this problem is obviously chaotic, but it is only over the 
past two decades that its connections^'*'^* to generic clean 
(i.e., without stochastic disorder) examples of quantum 
chaos^^ have been deepened. In both quantum chaotic 
and quantum disordered systems eigenstates are char- 
acterized solely by their energy, rather than by a set 
of quantum numbers. Since eigenstates and eigenvalues 
cannot be obtained analytically, it becomes useful to re- 
sort to some statistical treatment where one studies cor- 
relators averaged over large number of eigenstates instead 
of focusing on the properties of a single quantum state. 
Although the methods (and the language) of quantum 
chaos and quantum disorder have evolved independently, 
it was realized that Wigner-Dyson (WD) level statistics 
of RMT, a fingerprint^^'^° of quantum chaotic systems, 
is also applicable to disordered systems. ■^^'*'' However, 
the lack of any transport-related energy or time scale in 
RMT description of the delocalized phase signals that rel- 
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cvant time to traverse the sample diffusively, t-p = L"^ /D 
(T> = VF^Ii being the bare classical diffusion constant in 
3D), is set to zero in this framework. Therefore, it be- 
came clear that WD statistics can be applicable to disor- 
dered system spectra only for the energy separation scale 
much smaller than the Thouless energy''^ E^^ — fi/tj)- 

The physical origin of deviations from RMT in dis- 
ordered systems with 5 < oo is the finite time t-D re- 
quired for the particle to spread ergodically all over 
the sample, i.e., for the classical motion to explore the 
whole phase space. The statistical approaches to quan- 
tum disorder problems, like supersymmetric nonlinear 
(7- model (NLSM)^** which maps stochastic problem to a 
field theory without randomness, provide the justification 
of RMT in the ergodic {t ^ t-ri) regime of diffusive dy- 
namics. Furthermore, these techniques make it possible 
to study also the deviations from RMT for non-ergodic 
times or energy scales in weakly disordered {kpi ^ 1) 
conductors — both perturbative and nonperturbative cor- 
rections arc governed, within this framework, by the 
diffusion operator describing the dynamics of a corre- 
sponding classical system. The well-understood quan- 
tum chaotic properties of disordered system have made 
these systems a standard laboratory to test different 
approaches to generic quantum chaos^°'^^ (where well- 
defined averaging procedure over an ensemble is lack- 
ing^"'). Thus, limitations of RMT encountered in dis- 
ordered electron systems have put the study of devia- 
tions from the universality regime into the focus of both 
mesoscopic and quantum chaos communities, ^'^^ where 
"lessons from disordered metals" have greatly influ- 
enced the development of formalism for more arduous 
examples of quantum chaos. 

Recently, the equivalent program has been pursued for 
the eigenfunction statistics.'' The earliest prediction for 
the (universal) distribution of eigenfunction intensities, 
in Gaussian orthogonal ensemble (GOE) of random ma- 
trices (which are time-reversal and spin-rotation invari- 
ant), 



/pT(i) 
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:exp(-t/2). 



(3) 



is known as the Porther-Thomas^^ (PT) distribution. As- 
suming only that eigenfunctions are normalized but oth- 
erwise arbitrary, /pT(i) can be derived''** from the prob- 
ability that, e.g., component [which corresponds to 
^'(mo) at some point mo inside the sample] of an eigen- 
state of a N X N random matrix is equal to some value 
t/N 



fpT{t) = lim 



r{N/2) 



iv^oo N ^ T[{N -l)/2] 
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(4) 



in the limit N ^ 00 (which corresponds to an infinitely 
large lattice in our problem). This is essentially an exam- 
ple of the central-limit theorem, and it should describe 
completely the eigenstate statistics in the universality 
limit that is insensitive to any physical details of the 
system. Such limit in disordered electron systems re- 
quires infinite conductance g = Et^/A 00 since level 
spacing A (thermodynamic scale) sets the smallest en- 
ergy scale and t-p ^ <^ E^h 00 (in real systems 
i^Th is large only in small enough samples, such as quan- 
tum dots). The universality stems from the basis invari- 
ancc of RMT, i.e., the fact that eigenfunctions in RMT 
arc structureless with vE'a(r) and '^a{^') being uncorre- 
lated for |r — r'| > £, and fluctuating just as Gaussian 
random variables. However, random Hamiltonians of 
real disordered solids are tied to a real-space representa- 
tion, where matrix elements are spatially dependent and 
TBH (1) is a band diagonal matrix. Therefore, they do 
not satisfy statistical assumptions of the standard RMT 
ensembles since all elements of such random matrices are 
non-zero and spatially independent. Nevertheless, a rig- 
orous connection to the RMT eigenstate statistics [here 
just hcuristically established through the interpretation 
of Eq. (4)] is provided by Efetov's supersymmetric ap- 
proach^"'^^ (i.e., zero-dimensional limit of the NLSM). 
While WD statistics works well for the part of spectrum 
contained within the interval \E — E'\ <C -Exh, the distri- 
bution function f{t) in finite g systems (like the ones in 
Fig. 1) do not overlap with PT distribution in any interval 
of eigenfunction amplitudes. The redistribution of ampli- 
tude statistics, caused in part by the appearance of highly 
unlikely according to RMT prelocalized states, leads to 
three different regions of intensities t. The deviations are 
the strongest in the largc-t limit where f{t) can be orders 
of magnitude greater than PT distribution. This occurs 
also in quantum chaos^^ where localization due to scars 
is generally less pronounced than inhomogeneities of the 
prelocalized states in strong enough disorder. 

Small deviations of f{t) from the PT distribution are 
accounted by a WL-type correction^ (i.e., a quantity sen- 
sitive to the breaking of time-reversal symmetry; here we 
use the expression for GOE) 



fFuit) = fpT{t) 
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which is a regular function in the small parameter 1/g 
and can be derived by a perturbative treatment of the 
nonzero spatial modes of the NLSM. This result was ob- 
tained by Fyodorov and Mirhn* (FM) for t <C k^^^^ 
(where correction is smaller than the RMT term). As in 
the case of corrections to WD statistics, the deviations 
are parameterized by the properties of classical diffusive 
dynamics. Namely, k is defined in terms of the diffusion 
propagator (i.e., one-diffuson loop)^ 



Kdiflf = n(r, r) = 



gw' 



E 



exp(-Dq2r) 



(6) 
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which is the sum over the diffusion modes (diffusion prop- 
agator n(r, r') is the Green function of the diffusion equa- 
tion and can be expressed in terms of eigenvectors and 
eigenvalues of the diffusion operator — Pq^, in a rectangu- 
lar geometry studied here). The sum in Eq. (6) diverges 
linearly in 3D at largo momenta, thus requiring a cutoff 
at |q| ~ l/£to retain the validity of the diffusive approx- 
imation. We provide the ultraviolet regularization using 
exponential damping factor^^ which limits the sum to 
the diffusive regime Vq^ <^ 1/t {t = i/vp is the elastic 
mean free time) 
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Here the wave vectors q are quantized by the periodic 
boundary conditions used in all directions, = 2nnx/L, 
n2: = ±l,±2,..., and correspondingly for qy and q^. The 
argument y is expressed in terms of the semiclassical 
transport quantities 
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The sum (7) can be evaluated exactly by a simple nu- 
merical computation. On the other hand, its analyt- 
ical dependence on £/L is usually obtained after ap- 
proximating the discrete summation by an integral (a 
standard result in the literature is therefore quoted^ as 
S{y) ~ L/£). To avoid loosing the terms of the original 
sum, which depend on specific boundary conditions, it 
is necessary to evaluate the original discrete form. Fol- 
lowing Ref. 41 this can be done by using the function 



dSjy) 
dy 



-n^[2F{4y)f 



(9) 



Since F{y) is related to the complete elliptic integrals, 
for small values of argument y "C 1 they can be approxi- 
mated by the leading order to give 



ny) 



{T^y)~ 



1 



(10) 



which, upon integration in (9), leads to 



S{y) = ^+'-TT ln(y) - 3 tt^/^ ^ ^r^y + ap. (11) 

The integration constant ap can be fixed numerically by 
finding the limit, as y — > 0, of the difference between 
exactly calculated discrete sum S{y) and its analytical 
approximation (11) with removed constant ap. We find 
that this difference converges at small y < 10""' to ap ~ 
8.32. For the largest i in our study, the function (11) 
with this ap reproduces the numerically calculated sum 
over the diffusion modes to within 2%. Evaluation of 
S{y) in the standard way, by approximating it with an 



integral, gives only the leading term (which is the same 
for all boundary conditions) 



Si{y) 



ai, 



(12) 



where aj is an integration constant defined by the in- 
frared cutoff at small wave vectors |q| ^ L~^. 

The divergence of the sum over diffusion modes (7) in 
3D points out that more careful treatment is needed of 
the short-scale physics. Since k has the meaning of a 
time-integrated return probability for a diffusive parti- 
cle,^ it can be generalized to the ballistic case (various 
transport and thermodynamic phenomena encountered 
in disordered conductors are related to the classical re- 
turn probabilities for a diffusive particle, see Ref. 45). 
Using ballistic generalization^^ of the NLSM to go be- 
yond the diffusive approximation, it was shown ''■^^ that 
corresponding ballistic contribution to this return prob- 
ability, i.e., the probability for a particle to be scat- 
tered only once from an impurity and return back af- 
ter a time t <C r, has to be added to Kdis- Therefore, 
the total return probability is expected to be given by 
K = '^diff' + tbai- Such (scmiclassical) ballistic effects are 
non-universal, i.e., they can depend strongly on the mi- 
croscopic details of disorder (they are negligible in the 
case of smooth potential having a correlation length much 
bigger than'^'^^''* 1/kp). In fact, it was shown recently'^'^ 
that Kbai determining perturbative corrections in the An- 
derson model can be much greater than Kdis- In Sec. Ill 
we show explicitly that both contributions are needed to 
describe function of disorder strength, contrary to 

some previously reported results*^ for the 3D Anderson 
model where fitted k was related to ballistic effects only. 

The second point of intersection of computed distri- 
bution function with PT distribution can be considcired 
as the point at which a tail of f{t) starts to form. 
While fpuit) still approximately describes the "near tail" 
around these values of t, the slow decay of f(t) eventually 
requires a nonperturbative formula in 1/g. Such is pro- 
vided by the 3D case calculations*^ in the framework of 
standard diffusive NLSM, which give with an exponential 
accuracy 



/NLSM(i) ~ exp 



1 

4k 



t > K~ 



(13) 



and stems from the appearance of the prelocalized states. 
Since k ~ {kF£)~'^ (only in the leading order — the dis- 
crete sum in Eq. (7) also contains other terms), this re- 
sult written down to the leading order term of the cubic 
polynomial in ln[t/(fcpf)^]. 



ln/NLSM(i) ~ {kp£f In^ 



t 



{kpif 



(14) 



emphasizes that prefactor contains {kpi)"^ dependence on 
the disorder strength. However, the numerical constant 
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a 


i (a) 




W = 3.5 


21.5 


2.89 


0.00074 


W = 4.0 


17.4 


2.21 


0.0031 


W = 4.5 


14.3 


1.75 


0.0084 


W = 5.0 


11.9 


1.42 


0.018 


W = 5.5 


10.0 


1.17 


0.035 



TABLE I. Transport properties computed for our five en- 
sembles of impurity configurations characterized by different 
diagonal disorder strength W in the Anderson model on a 
simple cubic lattice 12^. 

is uncertain^^ because of being determined by the ul- 
traviolet q (in the field theoretical language), i.e., length 
scales ^ £ which are outside of the diffusive NLSM frame- 
work. Namely, the standard nonlinear cr-models^'*'^* are 
long-wavelength effective field theories for the diffusive 
modes whose saddle point is analyzed to get the cigcn- 
state statics. ""^^'^"'^^ The role of energy is played by the 
diffusion operator —W^, so that fundamental variables 
of the theory, Q{r) matrix fields, should vary spatially on 
the scale much larger than £. However, it was found^"'^^ 
that in 3D systems Q{r) vary rapidly on the length scale 

£, which then impedes the prospect of getting rigorous 
results and points out to a different physics determining 
the structure of cigcnfunctions in 3D, than is the case 
if low-dimensional systems. ^'^^ An attempt to overcome 
this limitations, by using ballistic NLSM which extends 
the semiclassical description to all momenta |q| < kp, 
leads to the same {kpiY prefactor dependence but offers 
7r/9-\/3 as the precise value of its constant piece. 

In between the perturbative {t <C k~^/'^) and the far- 
tail region {t > k^^) of the wave function amplitudes, 
NLSM analysis of f{t) predicts an intermediate range of 
amplitudes, described by^'^ 
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-t + 
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+ ■ 



(15) 



This has also has the form of the corrected PT distri- 
bution, as is the case of FM distribution. However, the 
correction term here is in the exponent, and therefore of a 
different type than the one in fpuit)- It should be large 
compared to unity, but small compared to the leading 
RMT term.^ 

An alternate route to account for the slow decay of f{t) 
at large amplitudes has been undertaken through a "di- 
rect optimal fluctuation method" (DOF) of Ref. 11. This 
approach suggests possible importance of ballistic non- 
semiclassical effects, which are missed in both diffusive 
and ballistic NLSM formulation, by analyzing the short- 
scale structure of solutions of the Schrodingcr equation 
(i.e., by analyzing the saddle point of the original problem 
of quantum particle in a random potential, rather than 



the saddle point of its effective field theory obtained by 
integrating out the disorder degrees of freedom). It gives 
the following result for the far tail asymptotics (where 
only the leading log-cube term of the full cubic polyno- 
mial of In t is computed explicitly) 



/DOF(i) ~ exp 

[-Cdof kF£ In^ t] 



(16) 



assuming Gaussian white-noise random potential, and es- 
timating Cdof — 3 • lO^'' for that type of disorder. The 
dependence on the prefactor on kF£ here is linear, which 
can substantially increase the probability to observe a 
rare event as compared to NLSM prediction. Nonethe- 
less, this result has also been interpreted heuristically 
within the NLSM picture,'' assuming that k might de- 
pend solely on the non-universal (semiclassical) ballistic 
^ l/kp£ contributions of the same kind as those encoun- 
tered in the region of small eigenfunction amplitudes. 



III. FITTING THE EIGENSTATE STATISTICS: 
PHYSICS BEHIND DEVIATIONS FROM 
RANDOM MATRIX THEORY 

In this section, we first describe the details of the com- 
putation of f{t), and also give elementary account of the 
quantum transport properties of conductors described by 
the Anderson model. This should serve as an overture 
for the subsequent detailed examination of the eigen- 
state statistics in the perturbative region (i.e., the main 
body of the distribution hmction of intensities), and non- 
perturbative region (where amplitude splashes generate 
large deviation from PT distribution), by attempting to 
fit the analytical forms introduced in Sec. II. We then 
analyze the confidence in such descriptions through rela- 
tive error of the corresponding fits. The extracted fitting 
parameters are interpreted by comparing them to the ex- 
pected ones in the analytical predictions, with the help 
of transport properties of our finite-size samples. 

We solve exactly the eigenproblem of the TBH (1) of a 
finite-size system inside a small energy window AE = 0.3 
positioned around E = Q. This interval picks 60 -70 
states in each of the 30000 conductors (modeled on the 
lattice 12"^) with different impurity configurations. The 
ensemble of disordered samples is characterized by the 
disorder-averaged conductance g. The overall number of 
collected states depends on disorder — as the disorder is 
increased the energy band broadens, meaning that some 
states start to appear with energy eigenvalues beyond the 
band edge Ei, = &t of the clean TBH while the average 
number of states in the band center decreases. These 
states are used to evaluate f{t) as a histogram of inten- 
sities at all points inside the sample {N = 12^). The two 
delta functions in Eq. (2) are approximated by box func- 
tions 5{x). The width /^E = 0.3 oi5{E) is such that p{E) 
is constant inside AE. The amplitudes of wave functions 
are sorted in the bins defined by 5{t — \^a(jc)\^V), where 
their width is constant on a logarithmic scale. 
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As emphasized in Sec. I, before embarking on the 
search for prelocahzed states, we first compute the 
transport properties of our samples. The exact zero- 
temperature conductance of the nanoscale finite-size 
sample is obtained from the Landauer-type formula'*^ 

g = Tr [t{Ep)t\EF)], (17) 

where transmission matrix t{Ep) is expressed in terms 
of the real-space Green fimctions^^ for the sample at- 
tached to two disorder-free semi-infinite leads. The de- 
tails of such calculations for the lattice model studied 
here arc given elsewhere. Here we just clarify the re- 
lationship of this conductance to the Thouless conduc- 
tance (/Th = 27ri?Th/A (expressed in terms of the spectral 
properties of a closed sample), which appears in standard 
analytical treatments of various disordered problems, in- 
cluding the eigenstate statistics.^ This mesoscopic com- 
putational technique opens the sample to the surround- 
ing ideally conducting medium, so that particles can leave 
and enter through the lead-sample interface. Thus, the 
discrete levels of an initially isolated sample are smeared, 
and spectrum of sample+leads=infinite system becomes 
continuous (which allows us to find the conductance at 
any Ep inside the band). However, the computed con- 
ductance, for not too small disorder^°'^^ or coupling to 
the leads of the same transverse width as the sample^°, is 
practically equal to the "intrinsic" conductance gxh- In 
practice this means that we attach sample to the leads 
of the same cross section and use the same hopping ma- 
trix element tmn throughout the system (i.e., in the leads, 
lead-sample coupling and in the sample), in order to min- 
imize any influence which leads can have on the conduc- 
tance.^^ 

The disorder-averaged transport quantities, charac- 
terizing the five ensembles of conductors studied here, 
are listed in Table I. For weak disorder W ^ 6 con- 
ductance g is dominated by the semiclassical effects. 
Thus, wc use the Bloch-Boltzmann formalism (applica- 
ble when^^ ^ 3> a), in Born approximation for the scat- 
tering on a single impurity, to get the elastic mean free 
path e.{EF = 0) ~ 35.4/VK2 shown in Table 1. Analytical 
treatments usually assume a simple spherical Fermi sur- 
face for which kp is just the radius of the sphere. Such 
quantity is not well-defined for a lattice system with non- 
spherical Fermi surface where kp is direction dependent 
(i.e., different average values can be obtained depend- 
ing on whether one averages the absolute value of kp 
or the root-mean-square of kp over the Fermi surface 
). Nevertheless, all different averaging procedures give 
similar values, and we use conventionally the one which 
would reproduce some transport formula, like Sharvin^^ 
classical point contact conductance G = Gqk^pL'^ /An, 
where such average values (here over the Fermi surface 
Ep = of a simple cubic lattice) naturally appear. This 
convention gives kp « 2. 8 /a, which should serve as a 
counterpart oikp appearing in theoretical simplifications 
assuming Fermi sphere. It is easy to check that these 



values of parameters, plugged into the Drude-Boltzmann 
formula g ~ (kpL^/An) (i/L), approximately reproduce 
the disorder-averaged Landauer two-probe conductances 
from Table I. 



A. Perturbative deviations from the RMT limit 

We commence the comparison between our results in 
Fig. 1 and functional forms from Sec. H by fitting FM 
distribution (5) to the data in the region of amplitudes 
where f{t) deviates only slightly from PT distribution. 
The plot of f{t)/fpT{t) - 1 (Fig. 3) offers a straight- 
forward way to check for the non-trivial feature of the 
correction term, such as whether the zeros Iq = 0.55 
and t'(^ = 5.45 of 3/2 - 3t + are exhibited by our 
data. We find that zeroes of the curves from our numer- 
ical simulation are slightly smaller (Fig. 3), the first one 
being even disorder strength dependent. Since naive (vi- 
sual) inspection of fits, especially on a logarithmic scale 
might lead to erroneous conclusions (like the range of 
amplitudes where some function fits the data), we es- 
tablish a quantitative criterion of the quality of a fit 
with some formula /fit(i) by looking at the relative error 
[f{t) — ffit(t)]/ f{t). This becomes especially important 
in assessing the fit of the far tail on the semilogarithmic 
scale (an example of such assessment is shown in Fig. 7). 
Moreover, this type of plot directly highlights intervals 
where some analytical formula can be considered to de- 
scribe the data. For example, this gives a quantitative 
insight into the boundaries of perturbative, intermediate 
and nonperturbative regions discussed in Sec. H. The FM 
distribution describes f{t) for small eigenfunction ampli- 
tudes completely in the weakly disordered samples, where 
[f{t) — fFM{t)]/ f{t) is less than 1% for the data and the 
fit in Fig. 3. It also fits the statistics in the conductors 
where dirty metal regime of transport is reached in the 
Anderson model (e.g., the resistivity of = 5.5 sam- 
ple would be ~ 500/if2cm for lattice spacing a = 3 A). 
However, it is demonstrated by Figs. 3 that increasing W 
toward the boundary W ~ 6 leads to larger relative dif- 
ference between the fitted FM distribution function and 
the numerical data, while also shrinking the interval of 
t where such comparison is still reasonable in the first 
place. For W ^ 6 unwarranted application of the Drudc- 
Boltzmann formula would give i < a (i.e., the transport 
in this regime is dominated by non-semiclassical and non- 
perturbative effects^^). Therefore, even though it might 
be possible to claim that FM form accounts for some 
portion of f{t) at stronger disorder,'*^ the direct com- 
parison of extracted parameter k to the calculated one 
becomes nonsensical. For example, W cannot be inter- 
preted as ~ in this regime. In fact, the departure 
of FM correction from our f{t) starts before the upper 
limit of disorder, determining the breakdown of Bloch- 
Boltzmann theory of semiclassical transport, is reached. 
Similar distinction between the disorder limitation for 
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FIG. 3. Weak localization correction k/2 [3/2 - 3t + t'^/2], 
through which /fm (t) accounts for the small deviations from 
the Porter- Thomas distribution /pT(t), fitted (dashed lines) 
to /(t)//pT — 1 in the perturbative region of amplitudes for 
statistics from Fig. 1. The first zero of all curves is W depen- 
dent and falls in the interval (0.46,0.51); the second one is at 
m 0.53. This should be compared to the two zeros t'o = 0.55 
and t[l = 5.45 of 3/2 - 3t + t^/2. This fit is expected to 
be valid from t < l/y^ to t = — we find the relative error 
[/(i)-/FM(i)]//(i) to be less than 1% all the way to t ~ 10"" 
(up to some noise in the data at the smallest investigated t). 



the diagrammatic results and somewhat higher hmit for 
the validity of Boltzmann equation is seen in some other 
cases^^ (despite the fact that Boltzmann theory is rigor- 
ously justified as a lowest order result of a perturbative 
expansion of e.g., the Kubo formula). 

The one-parameter fit of the FM result Eq. 5 allows us 
to extract the first relevant physical parameter from our 
data, K. The extracted values arc presented in Fig. 4. 
These values of k do not match to the calculated Kdis 
(Table I) determined by the universal properties of the 
diffusive dynamics. This is not surprising in the light of 
the discussion in Sec. 11. What might be surprising is 
that disorder-specific Kbai = k — Kdis can be few times 
greater than Kdis- This points out to the importance 
of the short-scale effects even in the perturbative region. 
Nevertheless, it appears that non-Gaussian features (non- 
zero higher order cumulants) of the uniformly distributed 
random potential of the Anderson model do not generate 
interesting terms beyond the second-order ones computed 
explicitly in the FM correction. 

For a system of fixed size. Figure 4 shows that k 
changes with the strength of disorder in a way ex- 
pected for the change of a sum of the two contribu- 
tions. The first term in the phenomenological formula 
for K = OMikpi)-^ + 0.22(A;f£)-i is different from the 
analytically computed Kdiff • However, the function fitting 
the data on the lower panel of Fig. 4 is only the leading 
order behavior (usually given when boundary conditions 
quantizing the diffusion modes are neglected*^) of the 
full Kdiir computed in Sec. II. Therefore, we proceed to 
fit formula (6) in the following form 
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FIG. 4. The extracted k from the fits of FM distribution 
(Fig. 3) to the portion of our numerical distribution f{t) 
contained within the interval of t where deviations from the 
Porter- Thomas distribution of RMT are small. The depen- 
dence of these values on the disorder strength (i.e., moan free 
path i) is explained as a sum of the diffusive contribution 
~ l/ikpef (or, more precisely, kdis expressed through the 
full sum (11), as shown in the upper panel) and the ballistic 
contribution ~ l/kpi- 



K = 



Ki2S{y) , K2 



(18) 



That such fit of k vs y = I/3(£/i)^ is successful is 
shown on the upper panel of Fig. 4, and also by com- 
paring K\/ y/y = 0.012L/£ factor in the first term to 
the expected l/g = {4TT/kl,L'^){L/e) = O.OII. The re- 
markably close numerical values confirm suggest that 
Kdiff + 0.23(fc_F£)~^ should include the lower order bound- 
ary condition dependent terms^^ in /tdifi- While it would 
be hard to investigate this extra terms from the scaling 
of numerically computed disorder-averaged conductance 
(because of "pollution" by the conductance fluctuations), 
here they parameterize the features of mesoscopic fluctu- 
ations themselves. 

The other important question concerns with the de- 
pendence of the two contributions on the sample size. 
Standard display of Kdiff and Kbai in the literature shows 
only the leading order terms of such parameters, like 
Eq. (14) or Eq. (16), which are L-independcnt.'''^^ The 
short-scale contributions to k are not sensitive to the 
sample boundaries, so that change of the fitted k in 
/FM(i) with L should be generated solely by a change 
of Atdiff(-^^)- Our evaluation of Kdis introduces several size 
dependent terms, cumulative effect of which is shown in 
the upper panel of Fig. 5. In the limit of large system size, 
this converges asymptotically to limL^oo Kdis{L) = 0.11. 
This offers a simple test of the accuracy of the regular- 
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FIG. 5. The region of small deviations of f{t) from the 
RMT statistics /pT(i) in the three sets of = 5 disor- 
dered conductors (lower panel), which differ by the lattice 
size: 12^, 16^ and 20^ (Fig. 2). Fitting WL correction 
k/2 [3/2 — 3t + t'^/2] to these curves gives k as a function 
of L, as shown in the inset of the upper panel. This should be 
compared to size-dependence of its diffusive contribution Kdis 
in Eq. (6), since ballistic contributions are size independent. 



ization scheme employed to obtain k(iis{L). In order to 
investigate this issue, we generate a set of conductors 
modeled on the lattices IG"* and 20'^ with W = 5 disor- 
der strength. The comparison between the correspond- 
ing distributions f{t) is shown in Fig. 2; using the same 
energy window AE = 0.3, this distributions are gener- 
ated from about ~ 150 or ~ 300 states, for each realiza- 
tion of disorder, which are picked in the band center of 
16"^ or 20'^ lattice, respectively. Even without too large 
statistics, a palpable deviation between the two distribu- 
tions is observed in the far tail, and also in the region of 
small t. In quantitative terms, analytical expression (6) 
gives Kdis{L = 12) = 0.019, KdisiL = 16) = 0.028, and 
Kdis{L = 20) = 0.036, while the fitted ones are k(L = 
12) = 0.089, k{L = 16) = 0.096 and k{L = 20) = 0.100. 
At lower disorder W = 3.5 we get somewhat better agree- 
ment between fitted and calculated change in k with the 
system size: k{L = 20) — k{L — 16) = 0.0036 versus 
Kdie{L = 20) - Kdis{L = 16) = 0.0027, in accord with 
the discussion above on the disorder strength boundaries 
for the validity of perturbation theory in l/kp^-- 



B. Far tail 

A common feature of (almost) all analytical results for 
the eigenstate statistics in 3D metallic samples is that 
some exp-log-cube formula is predicted to describe the 
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FIG. 6. Exp-log-cube fits, Eq. (19) and (20), versus at- 
tempted log-normal fit to the far tail of fit) for W — 5 ensem- 
ble of disordered conductors from Fig. 1. The numerical data 
for our three-dimensional systems clearly favor exp-log-cube 
asymptotics. Relative error of the fits is plotted in Fig. 7. 



large-t asymptotic behavior of f(t). Thus, the most gen- 
eral functional form along these lines would be an ex- 
ponent of a cubic polynomial of Int. Since lower order 
terms in the polynomial are rarely calculated, and fit- 
ting of any formula is more reliable when the number 
of free parameters is small, we choose to fit two different 
simple (employing only two parameters) exp-log-cube ex- 
pressions (postponing the physical interpretation of their 
parameters for a moment): 



(19) 



which is always the leading order term [and sometimes 
the only one amenable to explicit computation, like in 
Eq. (16)], and the NLSM-like result (13) 



/^,(t) = C;exp[- 



4k;elc 



In^(KELci)]- 



(20) 



Using the fitted k from Fig. 4, we first establish the 
boundaries of different intervals for f{t) which are dis- 
cussed in Sec. II (e.g., according to Eq. (13) the begin- 
ning of the far tails is approximately located at i > 1/k) . 
These serves as a consistent criterion, where portion of 
f{t) to be fitted with cxp-log-cube formula analytical en- 
larges with increasing W , thus avoiding spurious results 
when attempting to fit too large piece of the distribution 
function. Although log-normal asymptotics is expected 
in 2D systems^^"'^'^ (and has been confirmed numerically^^) 
it has also appeared as a candidate in 3D systems. How- 
ever, attempt to fit C^' cxp(— C2 In^ t) to our data fails 
completely, as shown in Figs. 6 and 7. 

To check quantitatively the level at which fitted exp- 
log-cube formulas match f{t) from the numerical sim- 
ulation, we perform our stringent test by plotting the 
relative error in Fig. 7. The possibility to fit both func- 
tions (13) and (20) in a more or less similar way is shown 
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FIG. 7. Relative error of the two different ways of fit- 
ting the exp-Iog-cube formula (as well as the illustration 
of inadequacy of the "log-normal fit" from Fig. 6, thin 
solid line) to the exact f{t) from Fig. 1: thick solid line 
is for /fit(t) = Cp cxp(— C3 In'^ t) and dashed line is for 
/4t(t) = Cpexp[-1/4K;BLC In^(KBLci)]- 



in Figs. 6 and 7. Thus, a portion of f{t) can be well- 
describocl by the oxp-log-cubc asymptotics, where the 
size of the interval of amplitudes where the relative error 
is small (because of the noise in the data points of f{t) 
being pushed to the largest values of t) increases with 
increasing disorder strength. Another conclusion which 
can be drawn from Fig. 7 is that a narrow intermediate 
region of amplitudes exists which is not covered by either 
the FM function or exp-log-cube asymptotics. However, 
it appears that it cannot be fitted at all by the intermedi- 
ate formula (15), which is too close to the PT distribution 
for this to work here. 

We now procecid by analyzing dependence of the ex- 
tracted parameters on the disorder strength kpi, as well 
as by looking at the consequence of interpreting kelc 
to have the same physical meaning as k obtained by fit- 
ting the FM formula to the perturbative region (Fig. 3). 
This interpretation is in the spirit of NLSM conclusions 
where both perturbative and nonperturbative corrections 
to RMT eigenlevel or eigenstate statistics are governed 
by the same semiclassical physics. However, their values 
turn out to be quite different along the respective por- 
tions of f{t). The most important finding in the far tail is 
that both prefactors C3 and 1/4kelc appear to increase 
linearly with kpt, as shown in Fig. 8. Since prefactors 
are the most conspicuous signature of the importance of 
underlying ballistic or diffusive effects, this would mean 
that ballistic effects completely dominate in the far tail 
(explaining e.g., why we find Kelc ^ k > Kdiff) and con- 
firming the conclusions of DOF analysis^^ where short- 



FIG. 8. Parameters of the two different exp-Iog-cube for- 
mulas [sec Eqs. (19) and (20)], fitted to the far tail of fit) in 
Fig. 1, as a function of disorder strength (measured by t). The 
success of the linear fitting of both C3 and 1/4kelc vs kri 
favors analytical predictions having first power of kpi as the 
prefactor of the leading log-cube order, thereby supporting 
DOF conclusions^^ which emphasize that short-scale effects 
govern the far tail in 3D systems. 



scale structure of the solutions of Schrodinger equation 
is pointed out to be solely responsible for the rare events 
at the largest possible intensities t. Nevertheless, a puz- 
zle remains: both prefactors depend on the sample size 
(which is not treated explicitly in the present analytical 
schemes^'^^), whereas it is plausible that ballistic effects 
are insensitive to L. For example, at = 5 we have 
to use C3 = 0.218 or kelc ~ 1.411 (compare with cor- 
responding values in Fig. 8) to fit the far tail of f{t) in 
Fig. 2 for the sample of size L = 20a. This suggests 
that a theory is needed where such size-dependent ef- 
fects are handled explicitly [with their magnitude cannot 
be explained by, e.g., the change in Kdiff(-f') as is approx- 
imately possible for the size dependence of parameters in 
the small-i region]. 



IV. CONCLUSIONS 

Our results confirm that statistical distribution of 
eigenfunction amplitudes in 3D mesoscopic disordered 
systems depends crucially on short-scale ballistic effects 
(and thereby on the details of a random potential). In 
fact, the exactly computed universal quantities character- 
izing the classical diffusion process generate much smaller 
contribution to the parameters needed to describe the ob- 
served distributions by the analj^ically predicted formu- 
las. This is revealed in the perturbative region of ampli- 
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tudes, where deviations of the eigenfunction amplitudes 
distribution function from the RMT hmit are small, and 
Fyodorov-Mirlin distribution (i.e., Porter-Thomas distri- 
bution corrected by the weak localization terms) captures 
their functional form for weak enough disorder (i.e., the 
upper limit of disorder strength should be smaller than 
the one determining the breakdown of other semiclassi- 
cal properties, like Drude-Boltzmann conductance) in the 
Anderson model of localization. In this interval of small 
eigenfunction amplitudes, the diffusive contribution con- 
tains boundary conditions dependent weak localization 
terms, which stem from evaluation of the discrete sum 
over the diffusion modes. In the region of large wave 
function amplitudes, deviation from the RMT appears 
in the form of far tail of the distributions function, and 
is governed by the prelocalized states formed in rare con- 
figurations of the random potential. The far tail, which 
we obtain after examining about 2-10^ exact eigenstates, 
cannot be accounted by the log-normal distribution, but 
is well-described by the cxp-log-cube formulas. However, 
substantially different parameters, used in NLSM formal- 
ism with clear physical interpretation, are needed for this 
when compared to the ones governing the small devia- 
tions. Nevertheless, their linear dependence on the dis- 
order strength (measured by kp£) is in accord with pref- 
actor of the DOF theory, where short-scale effects were 
pointed for the first time as the possible sole explana- 
tion of the far tail statistics (we cannot decipher whether 
these effects are of semiclassical or genuine quantum ori- 
gin). These findings appear to be quite different from 
the success of various scmiclassically based theories^'* in 
describing the spectral statistics in mesoscopic systems, 
which are dominated by the general properties of quan- 
tum coherent superpositions and diffusion. The problems 
with standard diffusive semiclassical description appear 
in the study of statistics of any quantity where correc- 
tions to RMT are expressed in terms of the sums which 
diverge above some specific space dimensionality^^ (such 
as d > 2 in the cigcnstate statistics problem). Our anal- 
ysis points out that regularization schemes, which ad hoc 
avoid divergences from the short-length scales below £ by 
introducing a cutoff at large momenta |q| ^ l/£, can lead 
to large discrepancies when compared to exact numerical 
results. 

The feature of the far tails in the Anderson model 
which does not fit into the picture of ballistic effects alone 
is the need for specific prefactors terms which can account 
for the observed size dependence of the statistics. One 
possible explanation for these discrepancies would be the 
lattice effects in our model, but real solids are lattice 
structures (here simplified by taking single orbital per 
site) and possibility of such effects would not be an arti- 
fact of the model (which has been a paradigmatic model 
of the localization theory since the seminal paper of An- 
derson^). The other possibility stems from the fact that, 
even though our samples are good metals with g > 10, 
they are not large enough to be in the vicinity 5 3> 1 of 
universality limit, which is the region treated by present 
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FIG. 9. Distribution function f{t') of the eigenfunction in- 
tensities t' = |^'(r)|^/|*(r)|typ (normalized by the typical 
value of intensity |^'(r)|fyp of a given eigenstate) in the An- 
derson modeled nanoscale conductors with W = 5 disorder 
strength on V = 12^ and V = 20^ lattices. This analy- 
sis is complementary to that shown in Fig. 2 for f{t) where 
t = |^'(r)|^/l/~^. The inset plots correction [analogous to the 
one in Fig. 3 for f{t)] to the FT distribution /pT(t') in the 
region of t' where such deviations are small. 



theories. Strictly speaking, our findings can be consid- 
ered as a demonstration of the structure of eigenfunctions 
in nanoscale conductors, and to contrast their structure 
to the predicted one in the limit g ^ 1, a, heuristic at- 
tempt to account for the small-size corrections might be 
useful. Thus, if we interpret t = |^'(r)py of the standard 
statistical analysis as the ratio of eigenfunction intensity 
to its typical value 1/V far away from the high ampli- 
tude splashes, then by the same token, in small sam- 
ples we can construct statistics of analogous quantity 
t' = |\E'(r)p/|\I'(r)|jyp, where eigenfunction intensity is 
normalized by the its typical value |^(r)|typ. The auspi- 
cious outcome of this procedure would be vanishing of the 
size dependence of the far tail statistics. It turns out that, 
for the same W = 5 examples (on 12"^ and 20^ cubic lat- 
tices) shown in Fig. 2, |^'(r)|jypy « 0.26 fluctuates only 
slightly from eigenstate-to-eigenstate or from sample-to- 
sample. This allows us to define the PT distribution of 
t', fMt') = [2nt'V\^{r%^]-'/^ exp(-t'y|*(r)|?yp/2). 
The correction term ,f(t')/fpT(t') — 1 in the region of 
small deviations from RMT now displays almost negli- 
gible size dependence (see inset in Fig. 9). Neverthe- 
less, Fig. 9 shows that size dependence of the far tail 
of f{t') persists, albeit with a smaller relative difference 
\fL=i2{t') - fL=2o{t')\/ fL=i2{t') than in the case of usual 
analysis dealing with f{t) in Fig. 2. 

It remains to be seen if present approaches can com- 
bine short-scale effects with those responsible for the 
size dependc;nc;c!. Namely, ballistic extension of the a- 
modcl,^^'^^ which attempts to overcome the shortcom- 
ings of the diffusion approximation and assumption of 
slow spatial variation of the NLSM fields, provides a pre- 
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cisc numerical coefficient of tlie log-cube prcfactor while 
leaving its dependence on disorder strength unchanged 
from (kpi)^ in the diffusive NLSM formalism this con- 
tradicts our findings. On the other hand, DOF tech- 
nique as formulated in Ref. 11 is not directly applica- 
ble to the conductors modeled by the standard Ander- 
son Hamiltonian.^^ Thus, the structure of eigenfunctions 
of the Schrodinger equation for a particle in a random 
potential (which is the basic, tantalizingly simply formu- 
lated one-particle quantum-mechanical problem) remains 
a problem to be elucidated further. 
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